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Abstract 



The recent measurement of = ^''"^ by the E821 Cohaboration at 
Brookhaven deviates from the quoted Standard Model (SM) central value 
prediction by 2.6a. The difference between SM theory and experiment may 
be easily accounted for in a variety of particle physics models employing weak 
scale supersymmetry (SUSY). Other supersymmetric models are distinctly 
disfavored. We evaluate for various supersymmetric models, including min- 
imal supergravity (mSUGRA), Yukawa unified 50(10) SUSY GUTs, mod- 
els with inverted mass hierarchies (IMH), models with non-universal gaugino 
masses, gauge mediated SUSY breaking models (GMSB), anomaly- mediated 
SUSY breaking models (AMSB) and models with gaugino mediated SUSY 
breaking (inoMSB). Models with Yukawa coupling unification or multi-TeV 
first and second generation scalars are disfavored by the measurement. 

PACS numbers: 14.80.Ly, 13.40.Em, 12.60.Jv 
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I. INTRODUCTION 



Recently, the Brookhaven E821 experiment has announced a new measurement of the 
anomalous magnetic moment of the muon [Q. The measured result deviates by 2.Qa from 
the central value of the quoted Standard Model (SM) prediction 0: a^(exp) — a^{SM) = 
43(16) X 10-1°, where also Ref. [^). The largest error on the SM calculation 

arises from the hadronic vacuum polarization loops. These loops are included via dispersion 
integrals involving the rate for e~^e~ — >■ hadrons, with the largest contribution coming from 
the region around the p(770) resonance.^] It is anticipated that improved measurements of 
these low energy cross section measurements will reduce the uncertainty in the SM predic- 
tion in the near future, without recourse to tau decay data. Furthermore, additional data 
already taken by the E821 collaboration should soon reduce the experimental uncertainty 
in by about a factor of two. If the difference between SM prediction and experimental 
measurement of is maintained and sharpened, then a clear signal for physics beyond the 
SM will be obtained. 

The error bars on are now small enough to be sensitive to electroweak loop effects, and 
other effects of the same order of magnitude. Weak scale supersymmetry is an especially well 
motivated extension of the SM in which contributions to from loops with supersymmetric 
(SUSY) particles naturally have a magnitude comparable to electroweak effects. Thus the 
Brookhaven experiment can potentially probe weak scale supersymmetry. Supersymmetric 
contributions to have been calculated previously [§-0 , including a number of very recent 
papers addressing the new E821 result |TD|-^]. The supersymmetric contributions to 
involve chargino-sneutrino loops and also neutralino-smuon loops. For tan/3 = — not too 
small, the supersymmetric contributions yield [0] 



where Mj (i = 1,2) is a gaugino mass, and Msusy is a characteristic sparticle mass scale. 
Then ^a^^^^ grows with tan /3, and for models with a positive gaugino mass, has the same 
sign as the superpotential Higgs mass term /x. 

In this paper, we present calculations of the contribution to from a variety of su- 
persymmetric models. Some of these models overlap with those recently examined in Ref. 
||TT| , p!3|Jl^ , p!6|JT9| , while others are new. In addition, we compare explicitly where possible 
with reach projections of Run 2 of the Fermilab Tevatron pp collider and the CERN Large 
Hadron Colhder (LHC). 

The supersymmetric models considered can be organized according to the assumed mech- 
anism for communication of supersymmetry breaking from the hidden sector to the visible 



^We note that other evaluations of the hadronic vacuum polarization contribution (primarily 
because of a larger uncertainty) can lead to a SM prediction in accord with experiment: see 
Ref. Q. The recent evaluation by Davier and Hocker |^ cited in Ref's ||l],^ uses the much more 
precise tau decay data to reduce the uncertainty in the e^e" — > vr+vr- cross section. This is 
a potential source of theoretical uncertainty since the tau decays purely via the weak isospin 1 
channel, while the photon also has an isospin zero component. 
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sector. First, we consider models with SUSY breaking communicated via gravitational in- 
teractions. These include the minimal supergravity model (mSUGRA), with universality of 
soft SUSY breaking terms at the Grand Unification Theory (GUT) scale. In addition, we 



also consider minimal SO{10) models with a high degree of Yukawa coupling unification p3 



models with an inverted scalar mass hierarchy (proposed to solve the SUSY flavor and CP 
problems), and models with non-universal gaugino masses. Next, we show results for models 
with gauge-mediated SUSY breaking (GMSB), for various numbers of messenger fields. We 
also examine models with anomaly-mediated SUSY breaking (AMSB) and finally, models 
with gaugino- mediated SUSY breaking (inoMSB). 

In general, for SUGRA and GMSB models, large CP violating phases may be present 
(such phases are suppressed in AMSB and inoMSB models). Phenomenologically, many 
of these phases must be small in order to satisfy constraints such as those coming from 
measurements of the neutron and electron electric dipole moments. However, it is possible 
that some CP violating phases are large, but still satisfy EDM bounds if certain amplitude 
cancellations take place. In this paper, we neglect the presence of CP violating phases on 
the g — 2 calculation; their effect on SUGRA models is noted in Ref. |^, and on general 



MSSM models in Ref. 10 



We end with some broad conclusions in the last section. If the apparent discrepancy 
between theory and experiment continues to hold, then SUSY GUT models with a high 
degree of Yukawa unification will be strongly constrained, as will be models which solve the 
SUSY flavor and CP problems via a decoupling of the first two generations of scalars. 



II. GRAVITY-MEDIATED SUSY BREAKING MODELS 

In this class of models, SUSY is assumed to be broken in a hidden sector, consisting 
of fields which do not interact with usual particles and their superparners via SM gauge 
or Yukawa type interactions. SUSY breaking is communicated to the visible sector via 
gravitational interactions. In general, there is no mechanism to suppress soft SUSY breaking 
terms which can lead to flavor changing (FC) and CP violating processes in conflict with 
experiment. A common (ad-hoc) solution is to assume universality of soft SUSY breaking 
terms at the unification scale, which suppresses the unwanted FC processes. 



A. Minimal supergravity model (mSUGRA) 



In this model it is assumed that, at the GUT scale, all scalars have a common mass 
mo, all gauginos have a common mass mi/2, and all trilinear soft SUSY breaking terms have 
a common value Aq. Electroweak symmetry is assumed to be broken radiatively (REWSB), 
leading to a model parameter set consisting of 

f^o, ^1/2, Aq, tan/?, sign{fi). (2.1) 

The superparticle masses and mixings can be calculated via renormalization group evolution 
between the scale of grand unification and the weak scale. We use the mSUGRA mass 
calculation embedded in the computer program ISAJET 7.51 for our results We have 
adapted the Aa^^^"^ calculation of Moroi for use with the ISAJET code. 
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Our results for the mSUGRA model are shown in Fig. We plot only parameter 
planes with /i > 0, since the opposite sign of /i almost always gives negative contributions to 
Aa^^^^. The solid shaded regions are excluded by i) a lack of appropriate REWSB, or ii) 
a lightest SUSY particle that is not the lightest neutralino (Zi), or in) by the experimental 
lower hmits from LEP2 that m^^ > 100 GeV, mg, > 100 GeV and m^^ > 76 GeV. The 
region below the thick solid contour has a lightest Higgs mass < 113.5 GeV, in apparent 
discord with recent results from LEP2 on searches for SM Higgs bosons. The dots 
shaded regions are found to have a value of Aa^^"^^ within 2a of the E821 result: e.g. 
11 X 10~^° < Aa^^'^^ < 75 X 10~^°. The regions below the dashed contours are accessible 
to SUSY searches via the isolated trilepton or signals at Run 2 of the Tevatron with 25 
fb~^ of integrated luminosity p5| , pp| . Finally, the region below the dot-dashed contour is 
accessible to SUSY searches at the CERN LHC pp collider, assuming 10 fb~^ of integrated 
luminosity ||30(] . 

In Fig. |I]a), we show the mo vs.mi/2 plane for tan/5 = 3, and for Aq = — 2mo. For larger 
values of Aq, we find the entire plane to be excluded| by LEP2 Higgs searches. At this 
low value of tan (3, the region preferred by the E821 result consists of a small region in the 
lower left corner. The preferred region lies entirely within the Tevatron Run 2 search region, 
but is unfortunately already disfavoured by LEP2 Higgs searches. In frame b), we show the 
same plane, but for tan/? = 10 and Aq = 0. In this case, the LEP2 Higgs bound cuts out 
a significant portion of the region preferred by Aa^^^^ , but a substantial region remains at 
low mo where ruh > 113.5 GeV. Much of this region is beyond the reach of Tevatron Run 
2 experiments, but all of the preferred region is well within the region accessible to LHC 
searches. Of particular interest is that models with large mo and low mi/2, i-e. in the focus 



point region |^ , are disfavored, and will be excluded if the disagreement between theory 
and experiment persists. Finally, in frame c), we show the same parameter space plane as 
in b), except now for tan/? = 35. The region preferred by Aa^^'^^ has expanded greatly, 
reflecting the nearly linear growth of Aa'^'^'^'^ with tan (3. Much of the region is beyond the 
LEP2 Higgs bound, but all of it is within the reach of the CERN LHC. In this case, the 
focus point region is allowed for large tan /3. 



In several papers, comparisons have been made between the calculated values of Aa 



SUSY 



and the b ^ s'y decay rate, and the neutralino relic density Q^^h"^. It is especially interesting 
to note that at large tan (3 and /i < 0, the mSUGRA model is doubly disfavored by both 



Aaj^ and by 6 — 37 ||3^. In addition, at large tan P, much of mSUGRA model parameter 



space at large values of mo and mi/2 is allowed by bounds on the neutralino relic density 
PB|^B[. This is due in large part to ZiZi ^ A,H ^ bb annihilation via very broad A 
and H poles in the s-channel. It is an important observation that the mSUGRA model can 
accommodate all these constraints, and simultaneously respect the bound from LEP2 Higgs 
searches. 



^Strictly speaking, the bound nih > 113.5 GeV applies to the SM Higgs boson and needs to be 
corrected for the SUSY case. However, over much of the parameter space h is close to the SM 
Higgs boson, and we will for simplicity use this as the limit in the rest of this paper. For this 
reason, as well as to allow for some uncertainty in the computation of m^, the reader should allow 
some latitude in the interpretation of this contour. 
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B. Yukawa unified SO{W) model 



Supersymmetric 5*0(10) grand unified models are especially attractive in that they can 
unify gauge couplings, Yukawa couplings and also matter particles within a single genera- 



tion. In light of evidence for neutrino mass from observation of atmospheric neutrinos [|37 
there has been heightened interest in this class of models SO (10) SUSY GUT models 
naturally accommodate light neutrinos via the see-saw mechanism. Using third generation 
fermion masses as inputs, Yukawa unification can be achieved if tan/3 ~ 50 and /i < 0. In 
this region of parameter space, however, REWSB is not possible assuming universal scalar 
masses at the GUT scale. Incorporation of D-teim scalar mass contributions, which are 
generically present when 5*0(10) breaks to a lower rank gauge group, allow REWSB to 



occur by imposing GUT scale boundary conditions where rriH^ < tuh^ In this case, 

the D-terms leave a characteristic imprint on the entire SUSY particle mass spectrum. The 
model parameter space consists of 

^^16! ^iQi Ml),mi/2, Aq, tan/3, and sign{fi). (2.2) 

Here, mig is the common soft SUSY breaking mass (renormalized at the GUT scale) of all 
matter scalars, while mio is the corresponding mass of the Higgs scalars. M^, parametrizes 
the magnitude of the 5*0(10) D-terms. Yukawa coupling unification restricts tan/3 ~ 50±4 
and /i < 0. The sparticle mass spectrum, relic density, 6 — >• 57 decay rate, and collider 
signals have been recently calculated in terms of these parameters in the second paper of 
Ref. m. 



One might expect these models to be heavily disfavored by the E821 result, since Yukawa 
unification occurs for /i < 0, while a positive value of Aa^^^^ occurs for /i > 0. We explicitly 
display these results in Fig. |^, for a) tan /3 = 47 and b) tan /3 = 50. We generate models with 
values of mi6 and mio over the range 0-3000 GeV, mi/2 : 0-1000 GeV, -3000 < Aq < 3000 
GeV and < 1000 GeV, and accept models that satisfy REWSB and phenomenological 
constraints, as well as having t — b — t Yukawa coupling unification at Mqut to 5%. In 
both frames, we see that Aaf^^^ is always less than zero. However, for tan/3 = 47, 



a 



substantial fraction of models have |Aa^^'^^| < 5, so that they lie within 3(T of the E821 
result. These models generally require mig > 1000 GeV, so are effectively entering the 
decoupling regime. In frame 6), for tan/? = 50, almost all models are excluded unless 
miQ > 2200 GeV. For the higher tan/? value, larger values of M2 and smaller values of 
are generated for Yukawa unified solutions; in some cases, neutralino loops dominate the 
amplitude for Aa^^'^^. Overall, for tan/3 = 50, there is decoupling at larger values of mie 
than in the tan/3 = 47 case. In these models, generally a reasonable value for ^-^Ji^ can 
be obtained Also, in these decoupling regimes, the rate for anomalous supersymmetric 



contributions to 6 57 decay rate also diminishes. However, such large values of mie may 
suffer from problems with naturalness. 

C. Inverted scalar mass hierarchy models 

A possible solution to the SUSY flavor and CP problems arises by decoupling the mat- 
ter scalars in the theory. This involves setting matter scalar masses towards the 10-100 
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TeV range, thereby suppressing all loop induced FCNCs and CP violating processes [HO 



Such models seemingly violate constraints from naturalness, but two ways out have been 
suggested. In the focus point region of mSUGRA models [^, large mo, small mi/2 and 
moderate to large tan (3 yields a small value for and possibly a small fine tuning of model 
parameters to achieve REWSB. An alternative is in models with an inverted scalar mass 
hierarchy (IMH), wherein first and second generations have multi-TeV masses, while third 
generation scalars have sub- TeV masses: this latter condition results in models satisfying 
"naturalness", since third generation scalar masses enter the fine-tuning calculation, while 
first and second generation masses are most severely constrained by flavor and CP violating 
processes. 



1. Radiatively driven IMH model 
An intriguing scenario has been proposed in Ref. PTl, wherein the IMH is generated ra- 



diatively by starting with multi-TeV masses for all scalars at the GUT scale. At large tan P, 
and for special soft SUSY breaking boundary conditions, third generation scalar masses 
are driven to weak scale values via renormalization group evolution, while first and sec- 
ond generation scalars remain at multi-TeV values, owing to their small Yukawa couplings. 
The proposed model incorporates t — b — t Yukawa coupling unification, and necessarily 
includes right neutrino superfields (as in 5*0(10)), and requires the special (S'O(IO) sym- 
metric) boundary conditions Aq = 2mfo = 4mfg for the renormalization group evolution of 
the soft SUSY breaking parameters. For consistency of Yukawa coupling unification with 
REWSB, 5*0(10) D-terms are again necessary Thus, the parameter space of the RIMH 
model consists of 

mi6, mi/2, tan/3, M|,, Mjy, and sign{fi). (2.3) 

The value of tan (3 is generally large to achieve Yukawa unification, and the parameter M^r 
is the scale of the superpotential singlet neutrino mass term: 10^ TeV < Mjv < Mgut, 
with Mat ~ 10^^ GeV favored by data from atmospheric neutrinos in the simplest seesaw 
model for neutrino masses. Only a modest IMH is possible. The mass spectra associated 
with these models, along with expectations for relic density, b ^ s'j decay rates and collider 



signals, have been discussed in Ref. |42,43 



Motivated by the considerations in Ref. P3|, we show in Fig. Oa) the masses of several 



16, 



sparticles, and in b) the associated value of Aa^ , for < 0, tan/5 = 50, Md = 0.2m 
"^1/2 = 0.25mi6 and Mn = 1 x 10^ GeV. These parameters give a reasonable IMH while 
maintaining a calculable sparticle mass spectrum. In a), the IMH is displayed by the mass 
gap between m^^ and m^^ and m^^. In b), we see that Aa^^"^^ is always negative for this 
case. However, for values of mig > 3000 GeV, the value of Aa^^^"^ is within 3a of the E821 
result, as the model moves into the decoupling regime. Such large values of mie appear to 
be beyond the reach of even the LHC, at least for an integrated luminosity of 10 fb^^ . 
We have also examined the situation regarding for /i > 0. In this discussed 



in Ref. P3| Yukawa couplings do not unify well, and the unification parameter R typically 
ranges between 1.6 and 1.9. Here, we have fixed = 0.2mie which tends to give the 
largest mass hierarchy H3] for the chosen tan/? = 50. The results of our analysis are shown 
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in Fig. ^ for a) Mn = 10^ GeV, and b) Mn = 10^^ GeV. The shaded regions are excluded 
because there is no REWSB, or one of the scalar mass parameters is tachyonic, or the lightest 
super symmetric particle (LSP) is not a neutralino. To the right of the solid line, the crunch 
factor S defined in Ref. exceeds 4. The dotted area shows the 2a region preferred by 
the result of E821 experiment. 

In Fig. ^, we show Aa^^^^ versus mig for the slice of Fig. with mi/2 = 0.22mi6 
for which the hierarchy tends to be large. We see that Aa^^^"^ is within 2a of the E821 
central value when rriiQ < 1700 GeV. SUSY signals should readily be observable throughout 
this range, which unfortunately is disfavoured by considerations of relic density and the 
branching ratio for the 6 — > 57 decay |[43||. 



2. GUT scale IMH model 

In this class of models, it is assumed the IMH already exists at the GUT scale. The 
model parameter space is 

mo(l), "^0(3), mi/2, Aq, tan/? and sign{fi), (2.4) 

where mo(l) is the common mass of first and second generation scalar s at Mqut, while 
mo(3) is the mass of all third generation and Higgs scalars at Mqut- In this case, two-loop 
contributions from first and second generation scalars to RG evolution of third generation 
scalars helps to drive the latter to small masses; to avoid tachyons, mi/2 must be chosen 
large enough. These models have been recently explored in Ref. [|^,^. In general, GUT 
scale IMH models can support much larger first and second generation scalar masses than 
RIMH models. This helps to further suppress FC and CP violating processes, but will also 
suppress anomalous contributions to a^. For instance, for the two case studies presented in 
Table 1 of Ref. we find Aa^^^^ = 1.6 x lO'^^ (case 1) and 1.8 x lO^^^ (case 2). Hence, 



these models should give values very close to the SM prediction. 

D. SU{5) models with non-universal gaugino masses 

In supergravity GUT models based on SU (5) gauge symmetry, gaugino masses arise from 
the Lagrangian term 

^ ^ (2.5) 

Mpi 

where the are the gaugino fields, and is the vacuum expectation value of the F 
term of a chiral superfield which transforms as the symmetric product of two adjoints under 

su{5y. 

(24x 24)syi„rnctric = 1 © 24 © 75 © 200. (2.6) 

Universal gaugino masses are obtained only if the superfield $ is a singlet of the GUT group. 
Higher dimensional $ representations yield non-universal gaugino masses at the GUT scale: 
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for the 24, Mi : Ma : M3 = -1 : -3 : 2, while for the 75, Mi : M2 : M3 = -5 : 3 : 1 and for 
the 200, Ml : M2 : M3 = 10 : 2 : 1. Thus, the model parameter space |46] is given by^ 



mo, M3, Aq, tan/3 and sign{jj,), (2.7) 

where the SU{2)l and t/(l)y gaugino masses can be calculated in terms of M3, the GUT 
scale SU{'i)c gaugino mass. 

These SU{5) models with non-universal gaugino masses will yield weak scale gaugino 
mass values very different from the mSUGRA prediction, and potentially also very different 
values of the muon anomalous magnetic moment. We have studied the SUSY contributions 
to within this framework for positive values of Ma/i. The 2a region favored by the 
E821 experiment is shown in the rriQ vs. M3 plane in Fig. ^ for tan/3 = 10 and Fig. |^ for 
tan/3 = 35. We fix Aq = in both figures. The frames a), b) and c) respectively show the 
cases where the superfield $ transforms as a 24, 75 and 200 dimensional representation of 
SU{5). The corresponding cases for the singlet $ are the mSUGRA cases in Fig. 06 and 
Fig. 0c. Again the shaded region in frames a) and c) is excluded by the same theoretical 
and experimental constraints discussed for the mSUGRA case. For the 75 case in frame b) 
the theory constraints together with > 100 GeV exclude the entire plane. However, in 
the 75 and 200 cases, the mass gap between the chargino and the LSP is very small, and 
the LEP constraint on the chargino mass may have to be reassessed. In these cases, the 
lighter neutralinos contain significant higgsino components, and LEP experiments may also 
be able to probe Z1Z2 production as discussed in Ref. [|^. In view of this, in frames b), we 



have chosen to include just the theory constraints in the dark shaded region. In the light 
shaded region (which covers the rest of the plane in Fig. |^), < 85 GeV. The various 
lines labelled by the sparticle type denote contours where mg^ = 250 GeV, = 250 GeV 
and rrig = 1,2 TeV. In frames c) the selectron is heavier than 250 GeV throughout the 
allowed region, while the gluino is always lighter than 1 TeV in the allowed range of frame 
b). Finally, the solid lines are contours where m/j = 113.5 GeV, except in frames b) where h 
is lighter than this throughout: in this case, the solid line labels rrih = 110 GeV. 

The area shaded by dots is the region favoured by the E821 experiment at the 2a level. 
We see that the 24 cases are qualitatively similar to the corresponding mSUGRA cases in 
Fig. 0, although the allowed region does not extend as far up in the GUT scale gluino mass. 
The bulk of this region would be accessible at a 500 GeV linear collider for tan/5 ~ 10, 
but even for the larger value of tan/5, a considerable portion of this area will be probed 
there. While simulations of LHC signals have not been performed within this framework, 
presumably this entire region will be probed by LHC experiments via the usual multi-jet plus 
multi-lepton signals. Multi-jet events with identified Z bosons and are a characteristic 
feature of this framework IHBI . 



Moving to the 75 case in frames b) , we see that the bulk (possibly all) of the parameter 
plane is already excluded by various constraints already described. Indeed if the region 
with charginos up to 85 GeV can be excluded either via chargino searches, or via neutralino 



^Here, we do not consider the possibility that an arbitrary linear combination of these irreducible 
representations is also possible. 
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searches, and h can definitely be determined to be larger than 110 GeV even after incor- 
porating effects of mixing amongst the CP even scalars, this case will be excluded. In any 
case, it should be possible to probe this entire region at future colliders, the near degeneracy 
of the Wi and Zi notwithstanding. 

Finally, in frames c), we see that the E821 favoured region extends over a large fraction 
of the plane for both values of tan/3. In the region to the right of the contour labelled 
Wi, rn^^ < 250 GeV, and reduces to below the experimental bound as we hit the excluded 
region. Chargino pair production will thus be accessible only over parts of the parameter 
space favoured by E821. We have checked that the lighter chargino is higgsino-like, or 
for relatively small values of Mg mixed, throughout the plane. Its detection may thus be 
complicated by the fact that the Wi — Z\ mass gap is small over much of this range. Again, 
although explicit LHC simulations have not been performed for this scenario, experiments 
at the LHC should be sensitive to the SUSY signal over a considerable fraction (if not all) 
of this parameter range. 

To get some flavour of the size of the various contributions to a^, we present in Table 1 a 
single case study for each of the models labeled by the dimensionality of the $ representation. 
The model parameters are also listed in the table. The corresponding sparticle spectrum 
is listed in Table 4 of Ref. |^. We recognize that this point may well be excluded by 
experimental constraints in some of the cases. Our purpose here is only to illustrate how 
the individual contributions depend on the dimensionality of $. In the table, we show the 
magnitude of the various SUSY loop contributions to Aa^^'^^, along with the total. In the 
case of 1 (universality), we see that the chargino loops dominate, while the Z2/ii [i = 1,2) 
loops are also large, but nearly cancel with each other. For this model, the Z2 is mainly 
wino-like. In the 24 model, at the weak scale Mi decreases by a factor of ~ 2 relative to 
the mSUGRA case, while M2 increases by about ~ 1.5. The loop contributions involving 
Wi decrease relative the mSUGRA case, but still dominate the total contribution. In the 
75 model, M2 increases by about a factor of 3 relative to mSUGRA, while Mi increases by 
a factor of about 5. The Z3 loops form the dominant neutralino contribution, and the total 
value of Aa^^^"^ is further suppressed. Finally, in the 200 case, M2 increases by a factor 
of about 2 at the weak scale relative to mSUGRA, while Mi increases by about 10. In this 
case, the Z4 is largely wino-like, and it gives the dominant contribution to Aa^^^"^ . 



III. GAUGE-MEDIATED SUSY BREAKING MODEL 

In recent years, there has been much interest in supersymmetric models where SUSY 



breaking is communicated via gauge interactions PSI. In these models, SUSY breaking oc- 
curs in a hidden sector, but SUSY breaking is communicated from the hidden sector to the 
observable sector via Standard Model gauge interactions of messenger particles, which are 
assumed to occur in complete vector representations of SU{5) with quantum numbers of 
SU{2) doublets of quarks and leptons. The messenger sector mass scale is characterized by 
M. The soft SUSY breaking masses for the SUSY partners of SM particles are thus propor- 
tional to the strength of their gauge interactions, so that squarks are heavier than sleptons, 
while the gaugino masses satisfy the usual "grand unification" mass relations, though for 
very different reasons. Within the minimal version of this framework, the couplings and 
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masses of the sparticles in the observable sector are determined (at the messenger scale M) 
by the parameter set, 



A, M, ng, tan/3, sign{fi), Cgrav (3.1) 



The parameter A sets the scale of sparticle masses and is the most important of these 
parameters. The model predictions for soft-SUSY breaking parameters at the scale M are 
evolved to the weak scale. The parameter Cgrav > 1 and enters only into the partial width 
for sparticle decays to the gravitino. 

In Fig. ^ we show the 2a region favored by the E821 data for the minimal GMSB model. 
Our plots are in the A vs. tan/? plane, for M = 3A, /i > and a) = 1, b) = 2, c) 
ns = 3 and d) = 2 but with /x = 0.75Mi. In all cases but d), the magnitude of /x is fixed 
by the REWSB constraint. The solid region is excluded either by the REWSB constraint, 
or by m^^ < 95 GeV, rrie^ < 100 GeV, m^^ < 76 GeV or < 100 GeV, as indicated by 
LEP2 searches. The solid contour denotes where rrih = 113.5 GeV. Other mass contours 
listed are = 250 GeV, = 250 GeV, m^^ = 250 GeV, and nig = 1 and 2 TeV 

(dot-dashed contours). The Wi, Ci and fi contours correspond to the approximate reach for 
SUSY particles of a Next Linear Collider operating at ^/s = 500 GeV. In frame a), we see 
that very little parameter space is favored for low tan /?, but considerable parameter space is 
favored for large tan f3. The stars correspond to the reach in A of the CERN LHC pp collider 
for specific GMSB model lines (assuming 10 fb~^ of integrated luminosity), with the tan/3 
values sampled in Ref. The large dots denote the corresponding reach for the Tevatron 



with 25 fb ^ of integrated luminosity 0]. In a), at least for low tan/3, even the Tevatron 



reach extends beyond the maximum value of A favored by the E821 result, while the reach 
of the LHC is far beyond. Similar plots are shown in b) and c) but for larger numbers of 
messenger fields. Frame d) is shown for the special case of a higgsino-like NLSP. For the 
model lines that have been analyzed |Q , we see that for the range of parameters favoured 



by the E821 data, SUSY signals might well be observable even at the luminosity upgrade of 
the Tevatron. 



IV. ANOMALY-MEDIATED SUSY BREAKING MODEL 



It has recently been recognized |^ that there exist loop contributions to sparticle masses 
originating in the super- Weyl anomaly, which is always present in supergravity models when 
SUSY is broken. In models without SM gauge singlet superfields that can acquire a Planck 
scale vev, the usual supergravity contribution to gaugino masses is suppressed by an addi- 
tional factor ^^^f^ relative to m| = M^^sy/Mp, or in higher dimensional models where the 
coupling between the observable and hidden sectors is strongly suppressed, these anomaly- 
mediated contribution can dominate. The gaugino masses turn out to be non-universal, and 
are found to be proportional to the respective gauge group /^-functions. Likewise, scalar 
masses and trilinear terms are given in terms of gauge group and Yukawa interaction beta 
functions. Slepton squared masses turn out to be negative (tachyonic). A common fix is 
to assume an additional contribution for all scalars. The parameter space of the model 
then consists of 

'^o? ''^3/25 tan/3 and sign{fi). (4.1) 
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In the minimal AMSB model (mAMSB), the Wi and Zi are both wino-like, and nearly mass 



degenerate, leading to a unique phenomenology |52[ . 

In Fig. 1^, we show plots of the mAMSB parameter space via the mo vs. 77^3/2 plane, 
for < 0, and for a) tan/5 = 3, b) tan/5 = 10 and c) tan/3 = 35. For AMSB models, 
/i < yields Aa^^'^^ > 0, since we take the gaugino masses Mi and M2 to be negative in 
ISAJET.0 Incidently, this sign of fi appears to be disfavoured by constraints on BR{b s-y) 
as obtained in Ref. [Q. The solid regions are excluded by lack of the correct pattern of 
REWSB, or when fi is the LSP, or by < 86 GeV (from LEP2 chargino searches in the 



mAMSB model ||5^). In frame a), only a tiny region at the lower left is within the 2a range 



of Aa^ as measured by E821. Moreover, in this entire plane rrih < 113.5 GeV. The solid 
contour denotes where rrih = 105 GeV: thus, even allowing for effects of mixing in the Higgs 
sector, much of this plane is presumably excluded. The dashed-dotted contour denotes the 
reach of the CERN LHC for mAMSB models with 10 fb~^ of integrated luminosity |^. In 
frame b), the solid contour denotes where rrih = 113.5 GeV, and the dotted shading denotes 
the 2a favored region of Aa^^"^^. Only a small region has Aa^^^"^ > 11 x 10"^'', and 
mfi > 113.5 GeV. We do not show any dashed-dotted contour as there was no computation 
of the LHC reach for this value of tan/?. In c), for tan/3 = 35, we find that a considerable 
region of parameter space is favored by the muon anomalous moment, while simultaneously 
having ruh > 113.5 GeV (solid contour). The reach of CERN LHC encompasses almost the 
entire favored region. 



V. GAUGINO-MEDIATED SUSY BREAKING MODEL 

A fourth class of models has been recently proposed, based on extra dimensions with 



branes, which provides a novel solution to the SUSY flavor and CP problems |£6[. In 
this framework, chiral supermultiplets of the observable sector reside on one brane whereas 
the SUSY breaking sector is confined to a different brane. Gravity and gauge superfields 
propagate in the bulk, and hence, directly couple to fields on both the branes. As a result 
of their direct coupling to the SUSY breaking sector, gauginos acquire a mass. The scalar 
components of the chiral supermultiplets, however, can acquire a SUSY breaking mass only 
via their interactions with gauginos (or gravity) which feel the effects of SUSY breaking: as 
a result, these masses are suppressed relative to gaugino masses, and may be neglected in 
the first approximation. The same is true for the A parameters. 

To gain a phenomenologically acceptable sparticle mass spectrum, it is necessary to 
include additional renormalization group running between the assumed compactification 
scale Mc and Mqut- At energies above Mqut, it is assumed that either S'0( 10) or SU{5) 
grand unification is valid. Thus, all scalars have masses mo = at Mc, but non-zero values 
at Mgut- The parameter space of the model consists of |H7 



mi/2, Mc, tan P, sign{^). (5.1) 



■^We have corrected an error in the sign of the A-parameters that was present in ISAJET v7.51. 
This does not affect the results in Ref. lls! 
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If one insists on a high degree of Yukawa couphng unification, then one is forced to require 
/i < 0. In this case, Aa^^*^^ will always be negative. Yukawa coupling unification also 
highly restricts the allowed values of tan (3. Other parameters may be necessary in addition 
to the above set depending on the specific model of grand unification assumed. 

In Fig. |10|, we assume the inoMSB boundary conditions to be valid at = 1 x 10^^ 
GeV, and that minimal SU{b) grand unification is valid between Mc and Mqut — 2 x 10^^ 



GeV. We use the SU{b) RGEs given in Ref. with Yukawa couplings ft = 0.519, and 
fi, = f^. = 0.277 at Mgut, which are characteristic of tan/3 = 35. Two additional SU{5) 
Yukawa couplings A = 1 and A' = 0.1 are assumed. Several representative sparticle masses 
are shown in frame a), along with the value of \fi\. Note that the lower bound on parameter 
space of mi/2 — 280 GeV occurs where m^^ = m^^. In frame b), we show the corresponding 
value of Aa^^"^^. As expected, it is always negative. However, for mi/2 > 1 TeV, it lies 
within the 3a bound from E821. In this case, the model is moving towards the region of 
unnaturalness. 



VI. SUMMARY AND CONCLUSIONS 

The recently reported measurement of by the Muon {g — 2) Experiment E821 has been 
interpreted as a harbinger of physics beyond the Standard Model. This will be so if in fact 
the quoted SM estimate of and its associated error are verified, and if the experimental 
measurement is maintained with a reduced error after the analysis of the year 2000 data. 
If the discrepancy between the measured and theoretical values of persists, then the 
result will act as a strong constraint on many forms of new physics, including weak scale 
super symmetric matter. 

We presented here the values of Aa^^"^^ expected in a variety of supersymmetric models. 
The quoted discrepancy |jl|] between theory and experiment can easily be accommodated in 
a variety of supersymmetric models. These include SUGRA and GMSB models with /i > 
and moderate to large tan (3. 

The SUGRA, GMSB and AMSB models at low tan (3 are disfavored. In addition, models 
with a negative value of 11M2 are disfavored; this includes models that incorporate a high 
degree of Yukawa coupling unification such as simple SU{5) or S'O(IO) SUSY GUT models, 
and inoMSB models. Finally, models that invoke a decoupling solution to the SUSY flavor 
and CP problems are disfavored. These include models with an inverted scalar mass hier- 
archy, and focus point models at intermediate {but not high) tan/3. Even so, like the SM, 
these decoupling models are generally allowed at the 3a level. 
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TABLES 



TABLE I. Loop contributions to Aa^^^^ in SU{5) models with non-universal gaugino 
masses. Each contribution must be multiplied by 10~^°. We adopt the parameter space point 
(mo, M3, Aq) = (100, 150, 0) GeV, with tan/3 = 5. For each case, /x > except model 24, for 
which /X < 0. 



loop 


1 


24 


75 


200 




44.5 


38.1 


11.6 


19.4 




-11.8 


-15.6 


-3.62 


-10.5 


Zifli 


4.58 


-1.57 


0.55 


-2.84 




21.2 


8.01 


0.71 


0.06 




-1.15 


-2.08 


-2.47 


5.49 


Z4/X1 


-0.25 


0.07 


-0.04 


48.5 


Zlfl2 


-7.56 


-1.32 


-1.31 


0.05 


Z2fl'2 


-21.8 


-9.74 


-0.69 


-0.37 


Zsfi2 


-1.49 


-0.77 


1.84 


-1.75 




6.14 


5.37 


0.96 


-33.3 


total 


32.4 


20.4 


7.53 


24.7 
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tan^=3, Ao=-2m(, tan^=10, A„=0 tan/?=35, A(,=0 




IDo (GeV) 

FIG. 1. A plot of mo vs. mi/2 parameter space in the mSUGRA model for |Lt > and a) 

Aq = — 2mo and tan/3 = 3, b) Aq = and tan/3 = 10 and c) j4o = and tan/3 = 35. The 2a 
region favored by the E821 measurement is shaded with dots. The region below the solid contour 
has m/j < 113.5 GeV. The region below the dashed contour is accessible to Tevatron searches with 
25 fb"^ of integrated luminosity, while the region below the dot-dashed contour is accessible via 
LHC sparticle searches with 10 fb ~^ of integrated luminosity. 
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FIG. 2. A plot of the Aa^'^^^ value for various Yukawa unified 50(10) models with a) 
tan/3 = 47, and fi < and b) tan /? = 50 and fi < 0. We require t — b — t Yukawa unification at 
Mqut at the 5% level. The pluses are valid solutions but have some sparticle or Higgs masses in 
conflict with LEP2 constraints. 
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FIG. 3. A plot of a) sparticle masses and b) Aa^^"^ versus mie for RIMH models with 
mi/2 = 0.25mi6, Md = 0.2mi6, tan/? = 50, < and Mn = 1 x 10^ GeV. 
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FIG. 4. The mie — mi/2 parameter plane of the RIMH model with Md = 0.2mi6, tan/3 = 50, 

II > 0, and a) Mat = 1 X 10^ GcV and b) Mn = 10^^ GcV.Thc shaded region is excluded by 
the theoretical constraints discussed in the text. To the right of the solid line, the crunch factor 
S > A. In the region shaded with dots, the SUSY contribution to is within 2a of the central 
value obtained by the E821 experiment. 
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FIG. 5. A plot of a) spaxticle masses and b) Aa^^^^ versus mie for RIMH models with 
mi/2 = 0.22mi6, Md = 0.2mi6, tan/? = 50, /x > and Mn = 1 x 10^ GeV. 
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FIG. 6. A plot of the mo vs. M3 parameter plane of the SU{5) model with non-universal 
gaugino masses discussed in Sec. IID for the case where the superfield $ transforms as a a) 24, b) 75, 
and c) 200 dimensional representation of SU (5). We have fixed ^0 = and tan /3 = 10. The shaded 
region in frames a) and c) is excluded by experimental and theoretical constraints discussed in the 
text. In frame b), however, the dark shaded region corresponds to just the theoretical constraints, 
while in the light shaded region < 85 GcV. In the white region around uiq ~ 500 GeV and 
= 125 GcV, the char gino is between 85 and 100 GeV. The various lines labelled by a particle 
type are contours of sparticle masses as discussed in the text. Finally, the 2a region favoured by 
the E821 experiment is shaded with dots. 
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FIG. 7. The same as Fig. except that tan/3 = 35. 
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FIG. 8. A plot of parameter space in the minimal GMSB model, for M = 3A, and > 0. 
In a), we take = 1, in b), = 2, in c) = 2> and in d), = 2 but with ^ = 0.75Mi. The 
2a region favored by E821 is shaded with dots. The reach in A of the CERN LHC for particular 
model lines with specific tan (3 values is indicated by asterisks. For instance, in frame a) the model 
line studied had tan(3 = 2 and the reach extending to A beyond 400 TeV. The large dots denote 
the corresponding reach of a luminosity upgrade of the Tevatron that yields 25 fb"^ of integrated 
luminosity. 
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FIG. 9. A plot of parameter space in the minimal AMSB model, for /x < and a) tan /? = 3, 
b) tan/3 = 10 and c) tan/3 = 35. In a), the solid contour denotes where nih = 105 GeV, while in b) 
and c) it denotes where nifi = 113.5 GcV. The 2a region favored by E821 is shaded with dots. The 
dot-dashed contours denote the boundary of the region that will be probed by the CERN LHC 
with 10 fb~^ of integrated luminosity. 
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FIG. 10. A plot of o) sparticle masses and 6) Aaj^^"^^ versus the minimal gaug- 

ino-mediated SUSY breaking model, for tan/3 = 35, and /x < 0. We assume SU{5) unification 
between Mc and Mqut- 
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